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Abstract 


An optical model description, based on multiple scattering theory , 
of longitudinal momentum loss in proton-nucleus and nucleus-nucleus 
collisions is presented. The crucial role of the imaginary component 
of the nucleon-nucleon transition matrix in accounting for longitudinal 
momentum transfer is demonstrated . Results obtained with this model 
are compared with Intranuclear Cascade (INC) calculations , as well as 
with predictions from Vlasov-Uehling-Uhlenbeck (VUU) and quantum 
molecular dynamics (QMD) simulations . Comparisons are also made with 
experimental data where available. These comparisons indicate that the 
present model is adequate to account for longitudinal momentum transfer 
in both proton-nucleus and nucleus-nucleus collisions over a wide range 
of energies. 


1. Introduction 

Propagation of a particle in the nuclear medium is 
characterized by its mean free path (refs. 1, 2, and 3). 
At incident energies of 50-150 MeV, the experimental 
mean free path for protons is known (refs. 4,5, and 6) 
to be approximately 5-6 fm. Theoretical calculations 
yield a value of similar magnitude only when Pauli 
blocking effects and the nonlocality of the optical 
potential are taken into account. These non-local 
effects can be satisfactorily accounted for only in 
the limit, where the imaginary part of the optical 
potential is small and all results can be expanded to 
first order in the imaginary potential (ref. 2). 

At high energies, however, the imaginary part of 
the optical potential is not small. In addition, the 
scattering cross section is extremely forward peaked. 
Thus, the concept of mean free path loses much 
of its significance as an indicator of how energy 
and momentum are exchanged. Alternate concepts, 
such as momentum degradation (or decay) length, 
have been introduced (ref. 7) in order to explore 
the inelasticities in nucleon-nucleon (NN) collisions 
at relativistic energies and their effects on nuclear 
stopping power. In nucleus- nucleus collisions, the 
first attempt to understand momentum degradation 
was based on a shock wave picture (ref. 7) where 
the nuclei were treated as colliding drops of “nuclear 
fluid” that formed a shock front where they met. The 
momentum degradation length A' in this approach is 
defined in terms of the conventional mean free path 

A = {Po<*Nn) 1 as 



where p 0 is the nuclear number density, &nn is the 
free nucleon- nucleon cross section, Pqm is the NN 
center-of-mass momentum, and is the mean lon- 
gitudinal momentum transfer per collision. These 
authors found that in general A' is an increasing func- 
tion of Pqm (fig- 1 of ref. 7). In particular, for bom- 
barding energies from the pion production threshold 
up to 600A MeV, A' remains approximately constant 
because of increased longitudinal momentum trans- 
fer and decreased elastic and quasi-elastic scatter- 
ings. Above 600A MeV, both Pqm and (^Q\\^ increase, 

but because they are nearly equal, A ; increases rather 
slowly. 

Further attempts to understand nuclear stop- 
ping power at high energies have replaced the shock 
wave picture with a two-fluid model (ref. 8), cou- 
pled only by a frictional drag. This model suc- 
cessfully demonstrated that high baryon density can 
be achieved during mutual interpenetration of pro- 
jectile and target nuclei without shock formation. 
The concept of momentum degradation length, how- 
ever, has survived the above transformation in think- 
ing about the nuclear stopping power. Anishetty 
et al. (ref. 9) examined the experimental data on pp 
(proton-proton) scattering and extracted information 
about proton stopping distance in high-energy colli- 
sions. Busza and Goldhaber (ref. 10) examined the 
proton-nucleus collision data and concluded that the 
mean rapidity loss suffered by a high-energy proton is 
2.4 ±0.2 units, which corresponds to a momen- 
tum loss of «4-7 GeV/c by the projectile nucleon 
in its rest frame. This is significantly higher than 
the momentum loss suffered by a proton in pp colli- 
sions (~0.75-1.5 GeV/c) (refs. 9, 11, and 12). The 
sequential scattering models (SSM) were introduced 
(refs. 13—17) in order to explain the significantly 



higher momentum loss in proton-nucleus collisions 
when compared with pp collisions. Csernai and 
Kapusta (ref. 17) defined the momentum degradation 
length in proton-nucleus collisions as 

A p = {Po&nnI) 1 (2) 

where /, called the inelasticity coefficient, plays es- 
sentially the same role as (p\\^ f PcM * n equation (1). 
This is not unexpected, since the mean longitudinal 
momentum transfer (p\\^ in AW collisions is inti- 
mately related to the degree of inelasticity in such 
collisions. A reasonable estimate (ref. 17) for A /; is 

a ”“G ^)^ 56fm < 3) 

where P z is the laboratory proton momentum. Note 
that not all versions of the SSM agree on this value 
of Ap, since they differ in their choices for the inelas- 
ticity coefficient I. The essence of these models can 
be illustrated by recalling that as an incident proton 
traverses a target, it loses some fraction of its momen- 
tum in each collision. Assuming a Poisson distribu- 
tion where the mean number of collisions is charac- 
terized by N(z) = Po&NN z (where z is the depth in 
femtometers) , one may write for P z 

P, - P„e- S £ (1 - If = (4) 

AT— 0 \ * / 

Hence equation (2) follows. 

In this regard, it is instructive to note that the 
rapidity loss distribution extracted (refs. 9-12) from 
pp collisions is Gaussian with a median rapidity loss 
shift of ^—0.7 and an average shift of «— 1 unit. 
In the rest frame of the incident proton, this trans- 
lates into a root-mcan-square momentum loss of 
~0.75 GeV/c (1.5 GeV/c if a variable proton to 
baryon ratio is taken into account). It has been noted 
(ref. 10) that the momentum loss suffered by the pro- 
ton in proton- nucleus collisions is roughly what it 
suffers in pp collisions times the number of inelastic 
collisions it is expected to undergo as it traverses the 
nuclear target. 

Although the above models describe rapidity loss 
of high-energy (>100 GeV) protons rather satisfac- 
torily, it is not obvious which model will apply in 
the low- and intermediate-energy (10 A MeV several 
A GeV) domains. It is in the low- and intermediate- 
energy regimes that we focus our present attention. 


The earliest evidence of energy and momentum 
loss in relativistic heavy ion (nucleus-nucleus) col- 
lisions in the few GeV domain comes from the 
observed “momentum downshifts” of projectile frag- 
ments in the pioneering experiments on relativis- 
tic heavy ion fragmentation using carbon and oxy- 
gen beams (refs. 18 and 19). In the projectile 
rest system, these momentum downshifts are small 
(typically tens of MeV/c). For charge-exchange 
channels, the measured momentum downshifts were 
larger (« 100 MeV/c). Gerbier et al. (ref. 20) re- 
ported a very large momentum downshift measure- 
ment of nearly 6 GeV/c for charge- pickup reactions 
of 900A MeV gold transforming into mercury in col- 
lisions with an aluminum target. This large down- 
shift, however, has not been reproduced in subse- 
quent experiments (ref. 21) where only small energy 
losses were measured («3-5 A MeV). Recently, Tull 
(ref. 22) measured the momentum losses of projec- 
tile fragments for a 1.65A-GeV argon beam collid- 
ing with both carbon and potassium chloride targets. 
Reference 23 contains further information on momen- 
tum loss in heavy ion experiments. Experimental lit- 
erature on stopping and compression in low-energy 
nucleus-nucleus collisions (<2A GeV) can be found 
in reference 24. Theoretical attempts to explain nu- 
clear stopping in this energy domain can also be 
found in the latter reference. 

The impetus for understanding momentum loss in 
proton-nucleus, as well as nucleus- nucleus, collisions 
within the multiple scattering theory framework of 
the present work was the desire to describe and pre- 
dict fragment momentum distributions in low and in- 
termediate energy (;=sl0A MeV to a few GeV) heavy 
ion collisions. The salient feature of this model is the 
association of the gradient of the imaginary compo- 
nent of the NN transition matrix (folded with nu- 
clear densities of the projectile-target system) with 
beam momentum loss (refs. 25 and 26). The real part 
of the transition matrix (folded with the nuclear den- 
sities) was previously shown to describe transverse 
momentum transfer (refs. 25 and 27). 

The outline of the paper is as follows. In section 2, 
the expression for dynamical momentum loss is in- 
troduced, and the analogy with the sequential scat- 
tering models is pointed out for the proton-nucleus 
case. In section 3, results of beam momentum loss 
are presented for 85A-MeV 12 C-induced reactions on 
targets ranging from 12 C through 197 Au. The im- 
pact parameter dependence of momentum loss for the 
reaction 12 C(84A MeV) 4- 12 C is also presented 
and compared with simulations and model calcu- 
lations, such as Vlasov-Uehling-Uhlenbeck (VUU), 
Boltzmann-Uehling-Uhlenbeck (BUU), Boltzmann- 
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Nordheim- Vlasov (BNV), and quantum molecular 
dynamics (QMD). Finally, in section 4 we conclude 
by summarizing the current status of model develop- 
ment and discuss future directions for experimental 
and theoretical research. 

2. Momentum Degradation 

In the optical model of momentum degrada- 
tion, one describes momentum loss in composite tar- 
gets within the multiple scattering theory framework 
(ref. 25) as a superposition of momentum losses by 
single NN collisions; that is, the latter is folded with 
the nuclear densities of the projectile- target system. 
For momentum loss in a single NN collision, the gra- 
dient of the imaginary part of the transition matrix 
is integrated from — oo to +z, where z is the depth 
(in fm), and the resulting expression is then folded 
with the nuclear density. Explicitly, 

Q\\ (M) = A p A t J (P^ppp J (PfrPr 

x j J V^plm i + J + — 

V) 

In the above expression, nuclear densities 
Pi (i = P,T) are normalized to unity, & ( i = P, T) are 
the nuclear internal coordinates, Ai are the mass 
numbers of the colliding nuclei, v = |u| is the relative 
velocity in the nucleon- nucleon center-of-mass frame 
with kinetic energy e = where p is the reduced 

mass. The gradient (V) is taken with respect to the 
projectile nucleus internal coordinates, and t is the 
complex two-nucleon transition amplitude given by 
equation (15) of reference 25. The bracketed quan- 
tity in equation (5) is the momentum loss in a single 
NN encounter, < 7 ||(r), with r = |r| the radial sepa- 
ration of the colliding pair. The imaginary part of 
the NN transition matrix is Im t. Note that Q\\(b, z) 

depends on both b and z> although asymptotically 
(z — > oo) it is only b (impact parameter) dependent. 

Now consider the transition matrix at high energy 

i (e, q) « (e) [a (e) + i] exp (- (e) g 2 ) 

( 6 ) 

where e is the NN kinetic energy in its center of 
mass and m is the nucleon rest mass. The root- 
mean-square momentum transfer, obtained from the 
real part of equation (6), is related to the imaginary 


component via an energy dependent ratio a(e), such 
that 



where a(e) is the ratio of real to the imaginary com- 
ponent of the NN forward scattering amplitude. The 
root-mean-square transverse momentum transfer is 

(<4al) = [5(e)] -1 = 0-3-0. 4 (GeV/c) 2 (7b) 

so that from equation (7a), one arrives at the follow- 
ing estimate for longitudinal momentum transfer: 

ydL) ~ °- 75 - 2 GeV / c ( 7c ) 

where an average of the experimental (absolute) val- 
ues for otpp{e) and a np (e), which range between 
as 0.2-0. 4 at high energy, was used. (The sub- 
scripts pp and np refer to proton-proton and neutron- 
proton, respectively.) The above estimate in equa- 
tion (7c) is consistent with experimental momentum 
degradation data in pp collisions (refs. 9 12). 

An important question is how the quantity Qy 
in equation (5) is related, in the high-energy limit, 
to the predictions of the high-energy models such as 
the sequential scattering model (SSM). Recall that 
while Q||(b, z) is both b and z dependent, geometry is 
eliminated in the latter models by explicit integration 
over the impact parameter so that only a tube con- 
taining N nucleons (N = 1, 2, 3, . . .) is encountered 
by the projectile nucleon. The number of nucle- 
ons encountered actually increases as the collisions 
become more central (i.e. , small impact parameter 
(b — > 0) collisions). Since an expression like Qj|(6, z) 
is explicitly density dependent, the present model 
offers some advantage over the high-energy models 
such as the SSM. 

Further insight can be gained by examining the 
Glauber model (ref. 28) expression for the aver- 
age number of nucleons encountered by an incident 
nucleon as it propagates through a nucleus: 

T ( b , z) « a NN f dz'pT (b, z') (8) 

J-OC V J 

This is shown in figure 1 for a proton incident on 
a lead target. We assume ~ 40 mb and use a 
Saxon- Woods density distribution for the lead nu- 
cleus with a half-density radius of 6.624 fm and a 
skin thickness of 1.73 fm. Also shown is the mean 
number of nucleons (encountered) in the SSM, that 
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z, fm 

Figure 1. Average number of nucleon's encountered by a pro- 
ton incident on a Pb nucleus as a function of impact pa- 
rameter b (fm), and depth along beam direction z . Also 
shown is mean number of nucleons in the SSM. 

is, N(z) = p 0 <TNN z » with Po = 0*17 fm" 3 , where z is 
the depth. Note that only for_srnalI impact param- 
eter collisions do T(b,z) } and N(z) agree with each 
other, as expected. 

3. Results 

In order to compare the predictions of this model 
with those from other models, as well as with ex- 
perimental results, we first focus on proton-nucleus 
collisions. The momentum loss suffered by rti inci- 
dent proton is plotted in figure 2 for a lead target 
at various depths z (fin) and at two different im- 
pact parameters as the projectile traverses the tar- 
get (the losses arc given in the projcctile-nucleus 
center-of-mass (CM) frame). The results are nor- 
malized to the asymptotic values of momentum loss, 
Q(6, oc), in each case. The maximum recoil momen- 
tum of the incident proton is shown in figure 3(a) as 
a function of target mass number, the latter rang- 
ing between 12 (for carbon) and 238 (for uranium). 
In figure 3(b), the mean number of inelastic colli- 
sions, V — ^T^inel/^reac i s plotted; ^mel and G reac 
are the proton- nuc leon i nel astic cross section and 
proton-nucleus reaction cross sections, respectively. 



z, fm 

Figure 2. Momentum loss as a function of depth by an inci- 
dent high-energy proton on a Pb target for two impact 
parameters. 

The value of <7j ne j at high energy is known to be 
«32 mb, whereas <r rea c is taken From compilations 
(ref. 29). Only an intermediate impact parameter 
collision is treated here; the actual magnitude of the 
impact parameter is target dependent. One perhaps 
might expect beam momentum loss to scale in a lin- 
ear way with v , the number of inelastic collisions the 
proton undergoes in the target. This is observed in 
our calculations as can be readily seen in figure 3(b). 

The number densities of target nuclei (Woods- 
Saxon type for A > 20 and harmonic well for A < 20) 
are taken from compilations (ref. 29). The appropri- 
ate equations to convert the quantities for the projec- 
tile in its rest system (PRS) into the center of mass 
(CM) or the laboratory (L) frame quantities and vice 
versa (CM to PRS, CM to L) are given in the ap- 
pendix. We are not aware of any set of experiments, 
other than the ones already cited, that measure recoil 
momenta in proton-nucleus collisions. This is not so 
in the heavy ion case, where such studies have a long 
history (ref. 30). Because of the paucity of data, we 
have decided to com par e our re sults with other theo- 
retical model predictions of longitudinal momentum 
transfer in proton- nucleus collisions. 
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(a) As a function of target mass number for a proton incident 
on targets ranging from 12 C through 238 U. 


(b) Versus mean number of inelastic collisions. 
Figure 3. Maximum recoil momentum of a high-energy proton (Ep = 1 GeV). 


One such model is the Intranuclear Cascade 
(INC) model and the code VEGAS (ref. 31). Results 
for the average forward momentum imparted to the 
struck nucleus as a function of bombarding energy 
are plotted in figure 4 for a 238 U target (also see 
fig. 2 of ref. 31). Results are plotted for the INC 
calculations and from this work. For the INC code, 
the average recoil momentum increases with bom- 
barding energy for the uranium target. It also in- 
creases with the mass number of the target (fig. 2 of 
ref. 31) in the INC. This is also true in the present 
work, as can be seen in figures 3(a) and 3(b). The 
energy dependence of the momentum transfer is sim- 
ilar in the present model, as shown in figure 4. Note 
that the present model yields a distribution of re- 
coil momenta of the target nucleus as a function 
of impact parameter, much as the INC does. Our 
calculation only yields an “average” value and, as 
currently formulated, is incapable of yielding stan- 
dard deviations around this mean value. Therefore, 
when a single value of recoil momentum is quoted, 
it is understood to be averaged over a range of im- 
pact parameters following reference 27, since the INC 
quotes an “average” value as well. Note that the av- 
erage forward momentum imparted to the struck nu- 


cleus is only a small fraction of the momentum of the 
incident proton, and this fraction decreases as inci- 
dent beam energy decreases. In figure 4 we plot the 
“average” value of forward (longitudinal) momentum 
(with and without pion production) obtained with 
our model. The value without pions ranges from 
89 MeV jc at an incident energy of E p = 400 MeV up 
to ^116 MeV/c at E p = 1800 MeV. From figure 4, 
the differences between predictions of our model and 
the INC are typically of the order of ~100 MeV/c 
or larger. Part of this may be due to the neglect 
of alpha and other complex fragment emissions in 
the INC code. Inclusion of these processes in the 
cascade would lower the average forward momentum 
imparted to the struck nucleus, thus reducing the dif- 
ferences between cascade predictions and the present 
work. The momentum “downshifts” (average paral- 
lel momentum transfer) associated with pion produc- 
tion are sizeable (refs. 26 and 32). It has been shown 
(ref. 32), for example, that the longitudinal momen- 
tum transfer (or the “downshift”) in pion produc- 
tion may range from ^370 MeV/c at threshold to 
«100 MeV/c at E p — l GeV in the projectile rest 
frame. These values, when converted to the lab- 
oratory frame, are of the order of the magnitude 
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Figure 4. Average forward component of momentum for 238 U 
target and predictions from INC calculations. 

of the discrepancy between the cascade predictions 
and those from this work (fig. 4). Further work 
is necessary to pin down additional sources of the 
discrepancy. 

For the nucleus- nucleus case, we first compare 
with the Intranuclear Cascade (INC) code of Yariv 
and Fraenkel (ref. 33), which is a direct general- 
ization of the VEGAS code (ref. 34). The code 
treats accurately the multiple collision processes in 
the nucleus but disregards completely the possible 
NN correlations (see refs. 33 and 34 for a detailed de- 
scription). The calculation reproduces the linear re- 
coil momentum distribution of the residual nucleus at 
the end of the cascade. The results arc shown in fig- 
ure 5 for 12 C + 12 C at E/A = 800 MeV. The histo- 
gram corresponds to the total recoil momentum. The 
differences may arise from the particular choice of 
physics input in the cascade code, for example, the 
lack of mean field or they may arise from the different 
choices for the nuclear densities in the two calcula- 
tions (harmonic well in the present work). Clearly, 
more work is needed in order to pinpoint the sources 
of these differences. 


Figure 5. Linear recoil momentum from present work and 
from INC calculations for I2 C + 12 C at E/A = 800 MeV. 


Next, we compare the predictions of the present 
work with those from the Vlasov-Uehling-Uhlenbeck 
(VUU) simulations by Aichelin and Stocker (ref. 35) 
and Aichelin (ref. 36) on stopping power for 84 A- MeV 
12 C-induccd reactions on targets ranging from 
12 C through 197 Au. At small impact parameters 
(6^1 fm), these authors found that for the 12 C tar- 
get, roughly 60 percent of the center-of-mass momen- 
tum was transferred to the midrapidity source, while 
projectile remnants continued with essentially beam 
velocity. This midrapidity source subsequently lost 
75 percent of its momentum because of NN colli- 
sions. For heavier targets, such as 58 Ni and 197 Au, 
roughly 66 and 80 percent of the maximum momen- 
tum transfer allowed by kinematics for each target 
was transferred. In addition, the projectile was com- 
pletely stopped in the target for Ni through Au tar- 
gets. In figure 6, we compare results from this work 
with those obtained by Aichelin and Stocker (ref. 35). 
The agreement is again reasonably good, especially 
for the heavier targets. 

The detailed impact parameter dependence of 
longitudinal momentum transfer is compared with 
the quantum molecular dynamics (QMD) calcula- 
tions of Aichelin (ref. 36) for 12 C + 12 C in figure 7 
for a bombarding energy E/A — 84 MeV. There is 
reasonable agreement. It can be seen from figure 7 
that the present model tends to predict longitudinal 
momentum transfers that differ by about 10 percent 
compared with VUU or QMD calculations. 
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Figure 6. Longitudinal momentum transfer per nucleon to 
target-like residue as a function of impact parameter b 
for 12 C-induced reactions on targets ranging from 12 C 
through 197 Au at E/A = 85 MeV. Triangles represent 
maximum momentum transfer possible from kinematics. 
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4. Concluding Remarks 

In this work, we attempted to answer the ques- 
tion, “How much of the incident momentum does 
a nucleon or nucleus lose as it collides with a nu- 
clear target?’ 1 This was done within the framework 
of a description of momentum loss based on multiple 
scattering theory. Using an experimental nucleon- 
nucleon (NN) transition matrix and the appropri- 
ate nuclear matter densities, we predict longitudinal 
momentum transfers to projectiles, targets, and their 
remnants. The beam energies considered range from 
hundreds of MeV to a few GeV per nucleon. These 
predictions are both impact parameter dependent 
and depth dependent. Results for proton stopping 
power are in reasonable agreement with the data and 
with predictions from the sequential scattering mod- 
els. For both proton-nucleus and nucleus-nucleus col- 
lisions, our predictions agree reasonably well with 
those from the Intranuclear Cascade (INC), Vlasov- 
Uehling-Uhlenbeck (VUU), and quantum molecular 
dynamics (QMD) model simulations. 

The capabilities of the present model are exten- 
sive. With only the imaginary part of the NN tran- 
sition matrix and nuclear number densities as inputs, 
the present model accurately predicts longitudinal 
momentum transfer in nuclear collisions. These pre- 
dictions depend on impact parameter 6, depth z, and 
energy. It is clear that just as the real part of the 
NN transition matrix largely determines transverse 
momentum transfer, the imaginary component pre- 
dominantly determines momentum loss suffered by 
the incident projectile on nuclear targets. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
November 3, 1993 


Figure 7. Longitudinal momentum per nucleon, carried by 
fragments in the reaction * 2 C + 12 C at E/A = 84 MeV 
as a function of impact parameter b (fm). Also shown 
are results from Aichelin (ref. 36) within BUU and QMD 
model frameworks. 


7 



Appendix 


Coordinate System Transformations 


In the rest system of the projectile, the center-of- 
mass velocity /3 qm and the Lorentz factor 7 cm are 
given as 


&CM = Pt/ (&T + M P ) ] 


^CM ^CM 


(Al) 


where Mp is the total mass of the projectile and Pp 
and Ep are the three momenta and energy of the 
target. 


In the present work, we calculate longitudinal 
momentum transfer in the projectile-target center- 
of-mass (CM) system. Let (Q z )> Pz > and P f z be 
the average longitudinal momentum transfer, and the 
initial and final momenta of the projectile in its rest 
system (PRS), such that 


(Qz) = P Z -P'= -p' z (A2) 


since P z = 0 by definition in this frame; ( Q z ) indi- 
cates an average value, since Q z (b ) is an impact pa- 
rameter dependent quantity, and we approximate it 
at one impact parameter b . The appropriate trans- 
formation from the center- of- mass frame (barred 
quantities) is then 

Pz = 7 cm (Pz - PcM M p) = -1 cm 0CmM p (A3) 
P'z = 7 cm (Pz - PCM E'p) = 7 cm (- (Qz) - (k.U E 'p) 


where E'p is the total energy of the projectile after 
the interaction in its rest system. A reasonable 
Ansatz for the latter is 


E'p — Mp + Q\/2Mp (A4) 

with Qp the transverse momentum transfer due to 
the interaction, which is also a function of impact 
parameter. 
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